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Galaxy interactions are a common phenomenon in clusters of galaxies. Especially major mergers are of particular impor- 
tance, because they can change the morphological type of galaxies. They have an impact on the mass function of galaxies 
and they trigger star formation - the main driver of the Galactic Matter Cycle. Therefore, we conducted a study of major 
mergers by means of a multi-phase ISM code. This code is based on a TREE-SPH-code combined with a sticky particle 
method allowing for star formation controlled by the properties of a multi-phase ISM. This is in contrast to the usually 
implemented Schmidt law depending mainly on the gas density. Previously, this code was used on isolated galaxies. Since 
our star formation recipe is not restricted to a special type of galaxy, it is interesting to apply it to interacting galaxies, too. 
Our study on major mergers includes a research of global properties of the interacting system, namely the star formation 
rate and the star formation efficiency, the evaporation and condensation rates, as well as the mass exchange of distinct 
components, namely stars, diffuse ISM, and clouds. Investigating these properties provides insight to interrelations be- 
tween various physical processes. The results indicate that the star formation efficiency as well as the evaporation and 
condensation rates are influenced by the interaction. 
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1 Introduction 

For some decades it has been known, that galaxy interac- 
tions play an important role in the evolution of galaxies 
changing their morphological type or triggering star forma- 
tion. Since it is possible to take re markable pictures of pe- 
culiar galaxies (see e. g. the question was raised 
if their features are a result of galaxy interactions. Peculiar 
features are e.g. bridges connecting galaxies or tails protrud- 
ing from them. The above mentioned assumption seemed 
likely, however, it was not clear until numerical simulations 
corroborated it. Th e first simulation s only considered gravi- 
tation al forces (e.g . Holmbe rdl94ll: Pleiderer & Siedentopf 
1961: lToomre & Toomrdl 1 972k iBarneslI 19881) . using differ- 
ent simplifications on the force calculation owing to the 
limited computing power at that time. Nevertheless, those 
simulations already demonstrated, that interactions between 
galaxies change their appearance. 

Simulations became more sophisticated with increasing 
computer power. Hence, it was feasible to address topics 
concerning the interstellar medium (ISM). Though the ISM 
only contributes approximately 10% of the mass of present- 
day galaxies, it is still an important ingredient: firstly, all 
stars are formed within the cold phase of the ISM, the mo- 
lecular clouds. Hence, it is necessary to consider this com- 
ponent in order to study the star formation rate during ga- 
laxy interactions. Secondly, the dynamics are strongly in- 
fluenced by hydrodynamics and especially energy dissipa- 
tion. Which properties of the ISM are considered depends 
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on the code used (see e.g. the review by lBarnes & Hernq uist 
Il992h . Our co de incorporates a 3 - phase ISM similar to the 
proposition by lMcKee & Qstrike il (ll977l) . The code was de- 
veloped in order to study Milky Way-like galaxies account- 
ing for several components, namely stars, dark matter, and 
the three phases of the ISM faarfst et al.ll2006l) . The warm 
and hot phases are treated by s moothed particle hydrody - 
namics implemented similar to iHernquist & Katzl dl989h . 
The molecular clouds are modelled by the sticky particle 
method (for details re fer to | Theis & Henslerll 19931) . A sim- 
ilar code was used by ISemelin & Combesl d2002l) . but they 
used a different sticky particle method and a different star 
formation prescription. 

Major mergers, that are interacting galaxies of about 
equal size, are studied by numerous authors. Often the re- 
sear ch's aim is to study the properties of the re mnant galaxy 
(e.g. lBarneslll988l: ISpringel & Hernquistll2005l) . or to look 
for substructures, such as tidal dwarf galaxies dDuc et al 



2004 iBournaud & Ducll2006l: IWetzstein et al.ll2007l) . Other 



goals are the r eproduction of observed major mergers (e.g. 



Karl et al.l l2008). or to study the overall properties of inter- 



acting galaxies, such as the star formati on rate (e.g. Mihos & 
Hernquist 1996: ldi Matteo et al. Il2008l) . While studying the 
star formation most often the S chmidt- Kennicutt relation 
dSchmidtlll959[ lKennicuttl[l998l) or some modification of it 
are used. The Schmidt- Kennicutt relation says that the star 
formation rate is proportional to the gas density, with the 
proportionality factor being a constant star formation effi- 
ciency. This work is distinct in that respect, since it does not 
use the Schmidt-Kennicutt relation, but uses a star forma- 
tion prescription following lElmegreen & Efremovl (11997b . 



®WILEY 

) imerScience® 



© 2009 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



790 



J. Weniger, Ch. Theis, & S. Harfst: A study of major mergers using a multi-phase ISM code 




Fig. 1 A scheme of all physical processes connecting the differ- 
ent components of a galaxy. In addition, the gas looses energy by 
radiative cooling and the molecular clouds by inelastic collisions. 



They found that the star formation efficiency is not fixed, 
but dependent on local properties of the ISM, i.e. the mass 
of the star forming cloud and the pressure in the ambient 
ISM. 

The main aim of this work is to study the evolution 
of the star formation during a major merger. In particular 
we are interested in understanding how a locally defined 
star formation efficiency affects the star formation rate. Our 
study also yields the mass exchange rates due to condensa- 
tion and evaporation. 



2 Interactions between different components 

The interaction participants are described by particles mim- 
icking distinct components of a galaxy, namely stars, dif- 
fuse gas, molecular clouds, and dark matter particles. Those 
components are connected by various physical processes 
(see Fig. [I]). The processes condensation, evaporation, and 
a drag force due to ram pressure allow for matter and mo- 
mentum exchange between the diffuse ISM and the clumpy 
molecular clouds. Star formation and related feedback of 
supernovae type II and planetary nebulae close the circuit 
of matter. Energy dissipation is caused by radiative cooling 
in case of the diffuse ISM and by inelastic collisions in case 
of molecular clouds. Since studying the evolution of star 
formation is the main aim of this work, the next section will 
describe its implementation in more detail. Fo r a thorough 
descr iption of all other processes we refer to lHarfst et al.1 
(200(1. 



2.1 Star formation 

The site of star formation is kn own to be giant molecular 
clouds (e.g. Lada & Lada|[2003). Since they are represented 
in our code by sticky particles, those particles have to be 
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Fig. 2 The star formation scheme according to lHarfst et all 
(2006). First the molecular cloud is inactive for a period n a (to). 
After that time span, a "star cluster" is formed within the molecu- 
lar cloud provided that all other star formation criteria are fulfilled 
(£i). The mass of the embedded star cluster is dependent on the 
local star formation efficiency, e. Finally, the cloud is fragmented 
due to energy input from supernovae type II and the ambient gas 
is heated (£2). 

converted to star particles. The implementation of it works 
as follows (Fig.0: 

- We assume that it takes ri a = 200 Myr till gas in a 
molecular cloud is able to form stars. During this time 
the molecular cloud is inactive. In other words, star for- 
mation is suppressed as long as the age of the molecu- 
lar cloud, £ now — ^cid is smaller than n a . ri a is the only 
free parameter in our star formation prescription and can 
be interpreted as a global star formation time scale. It 
is gauged by the star formation rate of the Milky Way. 
Since the mean star formation efficiency, e, in our sim- 
ulations is approximately 5% (see Fig. [4]) and the total 
mass in molecular clouds is approximately 4 • lO n M 
(see Table [B, it gives a mean star formation rate of ap- 
proximately SFR = rr 1 • e • M cki « 1 M© yr -1 . 

- Afterwards stars can be formed. They stay embedded in 
their parent molecular cloud. Their mass is determined 
by the star formation efficiency, e. A mass criterion and 
a pressure criterion prevent the formation of too small 
stellar particles: the mass of the molecular clouds, m c id, 
must exceed 2.5 • 1O 4 M and the pressure at the posi- 
tion of the molecular cloud, P gas , must be higher than a 
tenth of the ISM pressure at the position of the Sun in 
the Milky Way, i.e. P = 3 • 10 4 K cm -3 (Elmegreen 
& Efremov 1997). Note, that the star formation effi- 
ciency is not constant, but it depends on local proper- 
ties, namely the gas pressure and the mass of the parent 
molecular cloud. 

- After the most massive stars end their lives as super- 
novae, the molecular cloud fragments into four pieces 
of equal mass due to feedback by supernovae. 

3 Numerical setup 



Following lHarfst et al.l d2006|) the initial galaxies were con- 
structed by firstly ge nerating a disk/bulge/halo - system us- 
ing the method of iKuijken & Dubinski dl995 ). In the next 
step, one fifth of the stellar disk particles were transformed 
into molecular clouds and finally the diffuse ISM was added 
as an initially slowly rotating homogeneous sphere, which 
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Table 1 Mass of distinct components of one galaxy and 
corresponding number of particles 



Component 


Mass [1O 1A M ] 


NR. of particles 


Bulge- stars 


0.17 


10000 


Disk- stars 


0.29 


75098 


Clouds 


0.04 


36347 


Gas-particles 


0.02 


9017 


DM-particles 


1.94 


100000 


Total 


2.46 


230462 



will collapse in the firs t 200 Myr forming a warm d isk. The 
collisionless model of iKuijken & Dubinskil (I1995|) realises 
an equilibrium configuration. However, the equilibrium is 
affected by adding the ISM and the above mentioned pro- 
cesses. Therefore, we follow the system's evolution, until 
a quasi-equilibrium is established. The numerical integra- 
tion is done by means of a TREE-SPH code comb ined with 
the sticky particle method ( Theis & Henslerll 19931) . Gravita- 
tional forc es are determ i ned by the DEHNEN-Tree ( Dehnenl 
2002). See Ha rfst et al.1 (120061) for a more thorough descrip- 
tion. 

The relaxed galaxy has a mass of 2.46 • 10 11 M , thereof 
79% in dark matter particles. The remaining baryonic mat- 
ter consists of approximately 12% ISM particles - hence, 
molecular clouds and diffuse ISM - the rest being disk and 
bulge stars. The galaxy is realised by a total number of 
230462 particles. Note, that not only the mass, but also the 
number of clouds, SPH-particles, and disk stars as given 
in Table [T] will change throughout the simulations: this is 
caused by inelastic collisions, star formation, and feedback. 
The numbers of dark matter particles (10 5 ) and bulge stars 
(10 4 ), on the other hand, stay constant throughout the sim- 
ulations. Both galaxies are of equal mass and their rotation 
is prograde. The galactic planes coincide with the orbital 
plane. At t = the galaxies are placed at a separation of 
100 kpc. Their initial speed is chosen to match a parabolic 
orbit with a minimum separation of 20 kpc. The run was 
stopped at t = 3 Gyr. 



4 Results 

The code configuration allows to follow the evolution of dif- 
ferent components and their interrelations. To illustrate the 
time evolution, a series of snapshots are plotted in Fig. [3] 
exhibiting the surface density of the cold phase. The distri- 
bution of the cold phase resembles the distribution of the 
stars very well. This implies that dissipation does not affect 
the cold phase severely. In contrast, the diffuse ISM shows a 
slightly different distribution (see also Fig. [7]). The fist snap- 
shot of Fig. [3] (top, left) shows the initial configuration. The 
first passage at t = 400 Myr is pictured in the second panel 
(top, middle). Afterwards, a series of snapshots features on 
the one hand the expansion of the tails and on the other hand 
the development of a bridge, which vanishes gradually. The 
bottom row of Fig. [3] shows a snapshot shortly before the 



second passage which takes place at t = 1520 Myr, the be- 
ginning of the merging at t = 1650 Myr, and a quiescent 
phase after the merging (t = 2000 Myr). The beginning of 
the merging is defined by the last local minimum in the po- 
tential energy of the system, after which the centre of mass 
of the bulge components are no longer separating. In addi- 
tion, Fig. |4] shows the star formation rate, the star formation 
efficiency, the number of star formation processes per 10 
Myr, the evaporation and condensation rate and the evolu- 
tion of the masses of the different components. The vertical 
dashed lines in Fig.|4]indicate the first encounter, the second 
encounter, and the beginning of the merging. 

The star formation rate shows fluctuations due to the 
limited number of star formation processes during measure- 
ment intervals of 10 Myr. Still it is obvious that there is a 
maximum in the star formation rate at 1690 Myr, shortly 
after the merging starts. Another local maximum is given 
at 610 Myr, that is 210 Myr after the first passage. How- 
ever, Fig. |4] suggests an enhanced star formation already 
prior to the first passage. The period of enhanced star forma- 
tion located around the first passage is more than 400 Myr 
long. In this setup, the reason for enhanced star formation is 
not an enhanced number of star formation processes, but an 
enhanced star formation efficiency: the number of star for- 
mation processes does not correlate well with the star for- 
mation rate, but the star formation efficiency does, as one 
can easily infer from the left side of Fig. [4] Only the en- 
hanced star formation prior to the first encounter cannot be 
explained by an enhanced star formation efficiency. Here, 
it seems that the number of star formation processes is cru- 
cial. Although the mean star formation efficiency is only 
5%, the maximum star formation efficiency found within a 
time step can even exceed 50%. In general, Fig. [6] suggests, 
that the merging facilitates higher star formation efficien- 
cies. In fact, the maximum star formation efficiency is in the 
mean about 9% higher after t = 1650 Myr. Fig.[5]features a 
map of the star formation efficiency at t = 550 Myr, hence 
150 Myr after the first passage. The star formation efficiency 
is highest in the centre of the galaxies, but also sites of en- 
hanced star formation efficiency are recognised in the bridge 
and the tails. Though the star formation rate is affected by 
the interacti on, it is still not e xtraordinary high. Previous re- 
search (e.g. ICox et al.1l2004|) implies that the magnitude of 
the star formation rate depends on the initial conditions, es- 
pecially on the orbit. Furthermore the fraction of mass in 
molecular clouds might be too low to obtain a starburst of 
several tens up to several hundreds M yr _1 . 

Other important quantities are the evaporation and con- 
densation rates. The local maxima in the evaporation rate 
(at t = 630 Myr and at t = 1710 Myr) take place approx- 
imately 20 Myr after the local maxima in the star forma- 
tion rate. This time interval corresponds to the time delay 
between star formation and feedback due to energy input 
of supernovae. Note, however, that the time resolution is 
only 10 Myr. The condensation rate has its maximum at ex- 
actly the same time the merging starts. During this episode 
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Fig. 3 Snapshots of the interaction. Surface density maps of the clouds are shown (from left to right and top to bottom: t equals 0, 
400, 550, 750, 1000, 1250, 1500, 1650, 2000 Myr). The first passage takes place at t = 400 Myr, the second at t = 1520 Myr, and the 
merging starts at t = 1650 Myr. 
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Fig. 5 The star formation efficiency at t — 550 Myr. 
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Fig. 6 Evolution of the maximum star formation efficiency. 
The three vertical dashed lines indicate the first passage, the sec- 
ond passage, and the beginning of the merging. 



condensation outweighs evaporation. In general, however, 
evaporation exceeds condensation. The mean value for the 
former is 1 M yr _1 and for the latter O.8M yr _1 .A high 
condensation rate coinciding with the beginning of the mer- 
ging, has also been found in other major mergers performed 
by us so far. The understanding of this remarkable feature 
and its implications will be the topic of further study. 



Altogether the mass evolution of distinct components 
shows the following features (see middle and bottom dia- 
gram of Fig. |4]) 

- In general, the cloud mass is depleted by star formation 
while the stellar mass increases continuously. Further- 
more, the cloud mass decreases, because evaporation 
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Fig. 4 Evolution of star formation (left) and mass (right): the star formation rate (left, top), the evolution of the mean star formation 
efficiency (left, middle), and the number of star formation processes within a sampling period of 10 Myr (left, bottom). On the right 
side, the evolution of the evaporation and condensation rates (top), the evolution of mass of stars (middle), and the evolution of mass 
of the ISM (bottom) are plotted. The three vertical dashed lines indicate the first passage, the second passage, and the beginning of the 
merging. 



outweighs condensation. As a result only 3% of bary- 
onic matter are molecular clouds in the end of the sim- 
ulation. In comparison, the fraction of molecular clouds 
was initially 8%. 

- Due to star formation and related feedback during the 
first passage the stellar mass as well as the mass of the 
diffuse ISM rises more steeply. 

- As a consequence of the merging accompanied by en- 
hanced star formation, the stellar mass increases rapidly. 
In Fig. Hone can also distinguish a short period of time, 
within which the diffuse ISM is diminished by conden- 
sation, whereas the cloud mass increases. 

Fig. [7] shows a snapshot of the interaction at t = 550 
Myr - 150 Myr after the first passage. Different components 
are plotted separately in order to emphasise the multi-phase 
properties of our simulation. In the uppermost panel a sur- 
face density map of stars is shown. It is easy to pinpoint the 
two galaxies, a bridge connecting them, and tails protrud- 
ing from each of them. The cold molecular clouds (second 
panel from top) follow the distribution of the stars very well. 
However, the tail of the warm phase is considerably shorter 
as one can deduce from the third and forth panel which show 



a density map of the diffuse ISM and a temperature map, re- 
spectively. In the temperature map it is easy to identify two 
phases: a warm phase in the galactic disks, the bridge, and 
the tails, and a hot phase primarily located in the halo, but 
also in hot supernovae bubbles. The reason behind the short- 
ened tails is the initial setup of the diffuse gas component. 
Due to the initial collapse, the warm gaseous disk is initially 
much smaller than the stellar disk and grows only slowly in 
size. At the beginning of the interaction the warm gaseous 
disk is still smaller than the stellar disk. Since the utmost 
part of the tail evolves from the utmost part of the disk, the 
tail of the warm ISM is not as extended as the stellar tail or 
the tail of the cold phase. 

5 Conclusions & Outlook 

We have presented a study of an equal size maj or merger by 
means of a multi-phase ISM code described in Harfst et al. 
(2006). We have focused on the evolution of the star for- 
mation rate as well as on the evaporation and condensation 
rates. The multi-phase nature of the code allows for an ex- 
tensive analysis of mass transfer involving galaxy interac- 
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tions. Furthermore, the star formatio n prescription accord- 
ing to lElme green & Efremovl (I1997I) allows to follow the 
star formation efficiency during the interaction locally. The 
star formation rate as well as the star formation efficiency 
increases due to the first encounter and the merging of the 
nuclei. In the mean the star formation efficiency is only 5%, 
but at some locations, usually in the centre of the galaxies, 
but also in the bridge or the tails, the star formation effi- 
ciency can be much higher - up to 55%. Yet, a drawback of 
our models is an initial gaseous disk which is smaller than 
the stellar disk resulting in a lack of SPH particles at the 
tip of the tails. Therefore, an improved setup of the SPH 
component is needed in order to study substructures, such 
as tidal dwarf galaxies, within the tip of the tails. 
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Fig. 7 This figure shows from top to bottom a surface density 
map of stars, a surface density map of clouds, a density map of 
the diffuse ISM, and a temperature map of the diffuse ISM at time 
t — 550 Myr (i.e. 150 Myr after the first passage). 
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